Photoinduced charge and spin dynamics in strongly correlated electron systems 
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Motivated by photoinduced phase transition in manganese oxides, charge and spin dynamics 
induced by photoirradiation are examined. We calculate the transient optical absorption spectra of 
the extended double-exchange model by the density matrix renormalization group (DMRG) method. 
A charge-ordered insulating (COI) state becomes metallic just after photoirradiation, and the system 
tends to recover the initial COI state. The recovery is accompanied with remarkable suppression 
of an antiferromagnetic correlation in the COI state. The DMRG results are consistent with recent 
pump-probe spectroscopy data. 
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Competition among multiple phases is a key issue to 
understand electronic states in strongly correlated elec- 
tron systems. Nature of the phases is characterized by 
internal degrees of freedom of electrons, namely spin, 
charge, and orbital Low-lying fluctuations of these 
degrees of freedom grow up toward the phase bound- 
ary. Then, tiny amounts of external perturbations can 
dramatically change the electronic states. Prototypical 
examples are perovskite manganese oxides where a ferro- 
magnetic metallic (FM) phase competes with a charge- 
ordered insulating (COI) phase associated with antiferro- 
magnetic (AF) and orbital orders 0. In the COI phase, 
electric resistivity drastically decreases by magnetic field, 
pressure, impurity, and so on. 

Photoirradiation is known to be a powerful tool to in- 
duce such dramatic response. This is called photoinduced 
phase transition (PIPT) 0]. Ultrafast manipulation of 
femtosecond pulse laser is characteristic of the PIPT, and 
provides a striking difference between the PIPT and the 
other phase transitions. At early stages of the PIPT, 
only electronic part of the systems is strongly affected, 
and the spin, charge, and orbital degrees of freedom play 
major roles on relaxation dynamics. This is in contrast 
to the other phase transitions where structual deforma- 
tion sometimes occurs. Thus, we expect novel relaxation 
processes of the PIPT and electronic structures that can 
not be created in equillibrium conditions. 

The manganese oxides provide good playgrounds for 
elucidatingeffects of the multiple degrees of freedom on 
the PIPT [1H |E SI!, 03. In particular, the photoirra- 
diation into the COI phase strongly influences transport, 
magnetic, and optical properties. Remarkable decrease of 
electric resistivity has been observed in Pr .7Cao.3Mn03 
after the photoirradiation H H, @. Metallic behavior 
has been also observed in pump-probe reflection spec- 
troscopy measurements on Ndo.sCao.sMnOa 0]: spec- 
tral weight of optical conductivity above the COI gap 
is transfered into lower-energy region within 400 fs af- 
ter the photoirradiation. Magneto-optical Kerr spec- 
troscopy measurements performed on Ndo.sSro.sMnOa 
and Gdo.55Sro.45Mn03 suggest that a FM component is 
induced by the photoirradiation [j| Q . 

The purpose of this Letter is to examine microscopic 



processes of the PIPT in a system where the charge 
and spin degrees of freedom strongly couple with each 
other. Motivated by the spectroscopy measurements on 
the manganese oxides, we examine effects of the pho- 
toirradiation on the COI phase by calculating the tran- 
sient optical absorption spectra in the extended double- 
exchange (DE) model. Theoretical treatment of many- 
body interactions, dynamical response, and time depen- 
dence due to nonequilibrium states on an equal footing 
is challenging [if], [l2j]. In order to overcome this diffi- 
culty, we perform a density matrix renormalization group 
(DMRG) calculation [H M, M, QJ, [lj M, Gl ■ We find 
that the spectral weight above the COI gap is transferred 
into lower-energy region just after the photoirradiation, 
and then tends to go back to higher-energy region. This 
spectral-weight transfer is accompanied by remarkable 
suppression of AF correlations. In order to understand 
roles of localized spins on the transient spectra, the spec- 
tra of the extended Hubbard model are also examined. 
The DMRG results are consistent with the recent pump- 
probe spectroscopy data. 

The Hamiltonian of the extended DE model is defined 

by 

{ij),<r 1 (<j) 

-Jh ^2 ' C i,s(°)ss'Ci,s' + Js^S l - Sj, (1) 

i,s,s' (ij) 

where c\ a (cj )CT ) and Si represent a creation (annihila- 
tion) operator of a conduction electron with spin a at site 
i, and a localized spin operator, respectively. The elec- 
tron hopping, the on-site and nearest-neighbor Coulomb 
repulsions are taken to be t, U, and V , respectively. 
The AF exchange coupling betwwen the nearest-neighbor 
spins is taken to be J s (> 0). The conduction electrons 
interact with the localized spins by the Hund coupling 
Jh (> 0). We define 7ij = ri^j n-i.o- = c\ a Ci,<j, and 

Stot = J2i Si + J2i, s .s< k c l, s (^ ss ' Ci ' s ' witn tnc Pami ma_ 

trices a. The magnitude of the localized spin is assumed 
to be 1/2. From now on, the one-dimensional (ID) ver- 
sion of the model is taken into account. This simplifi- 
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cation is reasonable, since a ground-state phase diagram 
does not depend so much on dimensionality [20| . Ef- 
fects of the spin-charge separation on the photoinduced 
dynamics will be discussed later. 

Let us define the transient optical absorption spec- 
trum. First, pump light with energy luq comes into the 
system at time t$. At t\(= to + r > ^o)i a transient state 
is described by 
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where |0) represents the ground state with energy Eq, 

j = itJ2i,a(4,a c i+i,°- - H.c), Jq 1 corresponds to the 
dulation time of the pump light, and A is a normalization 
factor. Next, the transient state is probed by a photon 
with energy oj. The transient spectrum is given by 



a'(oj,T) = — hxx{ip(r) 
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where E\ = \ H \ip(r)) and 7 _1 corresponds to the 

dulation time of the probe photon. This form is valid 
for r ~7q~ + 7 _1 - In order to illustrate the effect of 
the photoirradiation on the optical spectra, a'(a>,r) is 
compared with the optical absorption spectrum a(u>) — 
(-l/7r)Im(0|j(^+^o--ff+i7) _1 j t |0). We also calculate 
O(r) = (^(t)| O \^{t)) for a static variable O. 

We apply the finte-system DMRG algorithm to cal- 
culations of a'(ui,T) and 0(t) [13]. An open bound- 
ary condition is taken for the odd number of lattice 
sites L in order to remove finite-size effects [2l|. Here, 
we take L — 13 and the electron number iV c i — 7. 
The density matrix of the system block is defined by 

Pi,j = J2k, a P a ^t*k^j,k^ where Y. a P a = 1 and k rep- 
resents a state in the environment block. A set of i\) a 

is taken Jo be {|0> ,jt | ) , |^(r)) , jt |^(r)) , (w + E x - 
H + i-f) 1 j' \iP(t))} for a calculation of c/(u>,t), and 
{|0},jt|0),|^(r)}}for O(t) 0, E ES HI M M ■ The 
truncation number is m = 256. A correction vector 
(u) + Ei — H + *7) -1 ^ |^( T )) is evaluated by the modi- 
fied conjugate-gradient method. The time-evolution op- 
erator e~ lHr is devided into n time slices e - lHr / n : an d 
the fourth-order Runge-Kutta method is applied for each 
slice. A value of n is determined so that tr/n = 0.01 cor- 
responding to E\t jn << 1. 

Let us first mention the ground state. The model 
parameters are taken to be U = lOt, V = At, and 
Jh = 2t being appropriate for the manganese oxides. 
The COI (FM) phase appears in a region that J s ^ 0A2t 

(J. < 0.12t) H3. A value of (0| |0> [= Stot(Stot+ 

1)] is minimum (maximum) in the COI (FM) phase. In 
the COI phase, spatial alternation of the charge density 
occurs, and the spin structure factor S(q) has a cusp at 
q ~ 7r. In the following calculation, we take J s = 0.2t, 



and the Hilbert space is restricted to $j?ot = 
Figure 1 shows a'(ui,r) for tr = 5 (= 
10 and 20, and a{u>) [23[. We take ujq 



0. 

To" 1 + 7- 1 ), 
= 2 At and 




FIG. 1: <x'{u, t) for various r values. The solid lines with and 
without the symbols represent a'(uj, r) and a(u>), respectively. 
The dashed line represents /3(u). 



7o = 7i = OAt [2J]. The optical absorption spectrum 
a(u>) has a finite charge gap being about 2t, indicating 
the COI phase. In contrast to a(ui), a'(ui,tT = 5) has a 
pronounced peak inside the gap at to ~ t. We have con- 
firmed that this peak shifts toward lower-energy region 
with increasing the system size L. Thus, the system just 
after the photoirradiation seems to be metallic. With in- 
creasing r, the peak shifts to higher-energy region and 
broadens. This shift suggests that the system tends to 
recover the initial COI state, and a part of the kinetic en- 
ergy of the conduction electrons flows into the localized 
spins. All of a'(u>, r) cross at uj ~ 2.8i, and for ui 2.8t, 
the spectral weight decreases with increasing t. We have 
confirmed that there is not remarkable difference between 
a'(uj,tT = 40) and a'{uj,tT = 20). 

In order to understand origin of the peak inside 
the COI gap, we calculate P(lo) defined by (3(u>) — 
(-l/7r)Im (0| i](uj + E - H + r^'V |0) with rj = rj - 

(0| rj |0) and 770 = ~tYli,a( c i,a c i+h<r + H.c). This repre- 
sents distribution of photoexcited states whose parity is 
different from that of a (u>). As shown in Fig. 1, the spec- 
tra a(uj) and j3{oS) are nearly degenerate for lo ^$ 3.5i, and 
f3(u) has a peak at u> = 3. It. The excitation energy of the 
peak inside the COI gap in a'(ui, tr = 5) is assigned to a 
difference between ojo (= 2.4t) and 3At, since the initial 
and final states of transition represented by a' (to, r) have 
the same parity as a(u>) and P(oj) have, respectively. 

In order to examine time evolution of the localized 
spins, we show E(t) — (ip{r)\ J2i Si ■ Si+i \^{t)) for var- 
ious J s values in Fig. 2. We find that E(t) increses 
with r and is saturated. Thus, the AF spin alignment 
in the ground state is disturbed by the photoirradiation. 
There are two characteristic time scales in E{t): (1) the 
dulation time of the pump photon, 7 ^ 1 , in which a differ- 
ence between E(0) and (0| J^i Si ■ Si+i |0) appears, and 
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FIG. 2: E(t) and (0| £\ Si ■ S i+ i |0) for various J s values. 
For all J 3 values, we take u)q = 2At for which a(u>o) has finite 
spectral weight. The arrows indicate tr values for which the 
transient absorption spectra are calculated in Fig. 1. 



(2) relaxation time of the localized spins that is about 
tr = 15 for J s = 0.2t. This value of the relaxation time 
is consistent with that estimated by the r-dependence of 
a'(uj,T) as mentioned in the previous paragraph. Thus, 
the charge and spin dynamics are highly correlated with 
each other. When J s decreases and the system ap- 
proaches the phase boundary between the COI and FM 
phases, the r-dependence of E(t) becomes remarkable. 
This means that the localized spins easily rearrange near 
the phase boundary. 




FIG. 3: a(ui) and a'(w,r) for the extended Hubbard model. 
The solid lines with and without the symbols represent 
a (u>, t) and a(u), respectively. The dotted line represents 
33a(a;)/13 for the extended DE model. 

Here, we make a comparison of a'(cj, r) between the 
extended DE and Hubbard models in order to illustrate 



roles of the localized spins on a'{u>, r). The ID extended 
Hubbard model is defined by 

H = -*£(4,«r*+l,*+ H -cO 

+U} j ni^n it i + V)^nin i+ i. (4) 

i i 

We take U = lOt and V — 6t in order to reproduce 
magnitude of the COI gap in the extended DE model. 
The other parameters are taken to be L = 33, N e \ = 17 
(9 electrons with up spin), and luq = 2 At. Figure 3 
shows ot'(u>, r) and a(u>) in the extended Hubbard model. 
We find that the COI gap collapses immediately after 
the photoirradiation. In contrast to the extended DE 
model, the r-dependence of a'(u>, r) is not remarkable, 
and the recovery of the initial COI phase is not observed 
for tr < 20. These differences of a'(u>,r) between the 
two models come from different characters of their spin 
correlations. In the extended Hubbard model, AF corre- 
lations are weak at quarter filling. On the other hand, in 
the extended DE model, the AF correlation between the 
localized spins is strong in the COI phase, and is remark- 
ably influenced by the photoirradiation. This change in 
the spin correlation causes the r-dependence of ot'(u>, r). 

We should stress the correlation of the spin and charge 
dynamics particulary in the extended DE model shown 
in Figs. 1 and 2. This correlation suggests that change 
in the charge states due to the photoirradiation strongly 
affects the spin states. This character of the transient 
states is in sharp contrast to the spin-charge separation 
for low-lying states in ID. This is because high-energy 
states with various charge and spin configurations are 
concerned in the early stages of the PIPT. 

Now, we discuss implications of the present DMRG 
results in light of the recent pump-probe spectroscopy 
measurements on the manganese oxides Hii. As al- 
ready shown in Fig. 1, a'(uj, r) in the extended DE model 
has the peak inside the COI gap just after the photoir- 
radiation. This peak shifts to higher-energy region with 
increasing r. These results reproduce the pump-probe 
reflection spectroscopy data [7|. It is also shown that 
this r-dependence of a'(u>, r) is accompanied by remark- 
able suppression of the spin structure factor S(q ~ ir) 
and slight enhancement of S(q ~ 0). Here, a finite value 
of magnetization does not appear, since the DMRG cal- 
culation has been performed on conditions S^ ot = (the 
minimum Stot) an( ^ \j> ^tot] = 0. These conditions may 
be reasonable in the manganese oxides where matrix el- 
ements of the spin-orbit coupling for the e ff -orbital van- 
ishes, but seem to be inconsistent with the Kerr rota- 
tion data [H, . For the electronic states in the man- 
ganese oxides, we suppose that short-range spiral-spin 
states or microscopic FM domains with keeping S* ot — 
are created just after the photoirradiation. Photocarri- 
ers can move inside these structures, while magnetiza- 
tion is not detected by the Kerr spectroscopy measure- 
ments whose spot size is of the order of 100 ^m x 100 
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/im. When the domain grows into macroscopic one, it 
is possible that volumes of up- and down-spin domains 
in the spot are different, and a finite Kerr rotation ap- 
pears. The above scenario is consistent with recent ex- 
periments on Ndo.5Sro.5Mn.O3 and Gdo.55Sro.4sMn03: a 
low-energy (high-energy) part of the reflectivity increases 
(decreases) immediately after the photoirradiation, while 
the Kerr rotation is delayed by the order of 1 ps [j| Q . 
The present DMRG data are expected to be applicable 
to the manganese oxides within a time scale before the 
Kerr rotation appears. 
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FIG. 4: Energy-surface diagrams. The filled and open circles 
represent Eo(St ot ) and Ei(S£ ot ), respectively. 

Finally, we examine stability of the macroscopic FM 
domain. Figure 4 shows an energy surface Ei(S* ot ) = 
(SU3HjUS^)/{S^\jj^\S^) as a function of S? ot 
that represents the energy of the photoexcited state just 
after the photoirradiation. The lowest-energy eigenstate 
of H with involved S£ ot is represented by \S* ot ). For com- 
parison, we plot E (S% ot ) = (Stotl-H" |Stot>- As shown 
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